Chaos in a one-dimensional compressible flow 
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We study the dynamics of a one-dimensional discrete flow with open boundaries - a series of moving 
point particles connected by ideal springs. These particles flow towards an inlet at constant velocity, 
pass into a region where they are free to move according to their nearest neighbor interactions, and 
then pass an outlet where they travel with a sinusoidally varying velocity. As the amplitude of the 
outlet oscillations is increased, we find that the resident time of particles in the chamber follows a 
bifurcating (Feigenbaum) route to chaos. This irregular dynamics may be related to the complex 
behavior of many particle discrete flows or is possibly a low-dimensional analogue of non-stationary 
flow in continuous systems. 

PACS numbers: 05.45.Ac 



Flows are a frequent topic of research among physicists 
- fluid flow traffic flow crowd movement [j], and 
granular flow Q are just a few examples. They are of 
particular interest because the particles that constitute 
the flow can exhibit complex behavior at certain flow 
parameters. There is a large body of work focused on 
understanding the cause of this motion and predicting 
the patterns and structures that these flows produce 111, 

nasi. 

Most flows consist of many mutually interacting de- 
grees of freedom and a complete description of the dy- 
namics is often impossible. It is not surprising then, that 
theories have historically focused on a statistical descrip- 
tion of complex flow [fj. Several studies, however, sug- 
gest that we can understand flows at a more fundamental 
level. Cellular automata fluid models are an example - 
for certain types of flows, the scale and particulars of col- 
lisions seem unimportant to the overall structure of the 
flowM, Q. The study of bifurcations in Taylor-Couette 
flow [9( suggest that the complicated motion in large scale 
flows can result from the interplay of a few chaotic de- 
grees of freedom. Several authors have successfully ex- 
tended these ideas to the general study of large coherent 
structures in fluid flow [lfjj. Finally, it has recently been 
shown that a one-dimensional series of nonlinear oscil- 
lators, when driven, can behave quite similar to larger 
scale turbulent flow (TlT | 

We seek to study complex open-boundary flow using a 
bottom up approach - by studying the simplest possible 
flows that exhibit complex motion. The low-dimensional 
model we present here is unique because particles inter- 
act with linear forces and the system has open bound- 
aries. There are many studies of low-dimensional, closed- 
boundary dissipative systems in the literature. Examples 
include driven, damped pendulums, the Lorenz equa- 
tions, and driven Frenkel-Kontorova models 13) 13 1 . Like 
these driven, dissipative systems, the system we present 
continually gains and depletes energy, but it also allows 
particles to pass across boundaries into and out of the 



flow. Research of complex, yet very low-dimensional, 
open-boundary systems is scarce in the literature. 

We present in this communication a simple one- 
dimensional flow with open boundaries that exhibits 
chaotic dynamics. The system consists of a line of point 
particles interacting with nearest neighbors according to 
a linear force law (Hooke's law). These particles travel 
towards an inlet at constant velocity, pass into a region 
where they are free to move according to their nearest 
neighbor interactions, and then pass an outlet where they 
are driven so that they have a sinusoidally varying veloc- 
ity. This outlet driving force continually supplies energy 
to the system and energy is dissipated when particles 
exit at any point away from their equilibrium position. 
As the amplitude of the outlet oscillations is increased, 
we find that the resident time of particles between the 
inlet and outlet follows a bifurcating route to chaos. We 
discuss the resulting dynamics of this system and suggest 
possible implications for larger dimensional flows. 
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FIG. 1: Diagram of the system. A series of point particles 
are connected by ideal springs and initially spaced S apart. 
Before reaching x = (the inlet) and after passing x = L (the 
outlet), each particle is constrained to the velocities shown. A 
particle moves according to its nearest neighbor interactions 
when between these points. 



The system consists of a one-dimensional chain of iden- 
tical point particles spaced a distance S apart and con- 
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nected by ideal springs. The particles travel in the re- 
direction and their positions are labeled Xi, where i is 
the index i = 1, 2, . . . , N. The number of particles N is 
chosen large enough so that the flow is sustained through- 
out our simulations. Particles undergo different dynam- 
ics as they pass certain points along the flow. Fig. Q] is a 
schematic of the system. 

The inlet is located at x = 0, the outlet at x = L, 
and the region in between is referred to as the chamber. 
The time when particle i reaches x = is labeled t\ n 
and is calculated t\ n = iS/v. The time when particle i 
reaches x — L is labeled t l out and is determined implicitly 
from the equation Xi(t l out ) — L. The velocity of a par- 
ticle before and after these times is constrained as given 
in Eqs. Ijlall bp . Particle velocities are not constrained 
between t]„ and t\„ t . 



±i = v t < t\ n , (la) 



x i =v + AcoB(u(t-t % out )) t>e out . (lb) 
Here A and ui are the amplitude and frequency of velocity 



oscillations after reaching the outlet. 

The initial spacing between particles is S, a value cho- 
sen such that the time between particles entering the 
chamber is much greater than the average resident time 
of individual particles within the chamber. This ensures 
that only one particle is unconstrained at any moment in 
time. 

The following are the complete equations of motion for 
the particles. 
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Xi = vt — iS t < t 
k(xi-i + Xi+x) 



W + Ac0B(w(t-O) 
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(2a) 

„ut, (2b) 
(2c) 



All particles have the same mass, m, and the linear 
restoring force, k, is the same for all springs. The solution 
for the position of the ith particle for times t\ n < t < 
t % out is given in Eq. ([3|). It is used to solve Xi(f out ) = L 
implicitly for t l out using Newton's method. 
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Xi (t) = Bi cos(V2ai) + B 2 sin(V2ai) 
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where, 
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In all of the following results, we use the following pa- 
rameters: S = 100, L = 2, v = 10, a = .06, u = 13.2, 
A = to 40. These parameter values ensure that, at 
most, only one particle is located in the chamber at any 
moment in time. 

For each particle we calculate the resident time within 
the chamber. This quantity is defined, 



I ... 



(5) 



t\ n is determined from Eq. [Ta| and t l out is determined by 
numerically solving Eq. [3] for Xi(t l out ) = L. 

Fig. [3 shows several return maps for the resident 
time (these map t l res to t 1 ^}). In general = 
f(t l res7 t l r ~s, ...), but for this system, 



4 + 1 _ 
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(6) 



i.e., the return map is one-dimensional. This results from 
constraining particle velocities before x = and after x = 
L - reducing the system to only one degree of freedom. 
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FIG. 2: Return map for t\ ea for several values of the am- 
plitude, A. Parameters for this simulation were: S = 100, 
L = 2, v = 10, a = .06, u = 13.2. 
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FIG. 3: a) Bifurcation diagram of t\. Ba plotted for values of A from 25 to 40. b) The Lyapunov exponent, A, as a function of the 
amplitude of outlet oscillations, A, also from 25 to 40. Parameters for this simulation were: S = 100, L = 2, v — 10, a = .06, 
Lu — 13.2. The first three bifurcation points are located at A = 29.0, A — 32.1, and A = 33.1; and A > first at A = 33.4 



As A is increased, a peak develops in the return map 
and the fixed point eventually becomes unstable. The 
one-hump map that develops is likely to exhibit the fea- 
tures of many other unimodal return maps - specifically 
a bifurcation route to chaos Q. As Fig. [3^, shows, this is 
indeed the case. The system is iterated onto the attrac- 
tor and the next 100 values for t l res are plotted for values 
of the amplitude, A, from 25 to 40. The resident time 
bifurcates several times and eventually becomes chaotic 
before settling back to a period three dynamics. The 
first three bifurcation points are located at A = 29.0, 
A = 32.1, and A = 33.1. 

In Fig. [3j3 we plot the Lyapunov exponent A for the 
resident time as a function of the amplitude A. The 
Lyapunov exponent is a measure of the separation of in- 
fmitesimally close trajectories and in this case is calcu- 
lated numerically from the following equation, 



1 n—1 

km -Vlog|/'(4 es )|. 



(7) 



When A > trajectories exponentially diverge, which 
produces chaos when the trajectories remain bounded. 
The system becomes chaotic at A = 33.4 where A first 
turns positive. 

The system we present above is quite different than 
other one-dimensional particle models in the literature 
111 1 1 31 ] . Instead of using nonlinear interactions between 



particles, the particles in our system interact with lin- 
ear forces and constraints are applied abruptly at the 
boundaries. This shows that complex motion can arise 
in a flow at the boundary between simple constrained 
motions without the need for nonlinear interactions be- 
tween particles. Many large scale flows contain regions 



where the dynamics are tightly constrained to regular 
motion, with complex motion occuring at the boundaries 
between these regions. Simple models such as the one we 
have presented can provide insight into how this behavior 
develops. 

Summarizing, we have presented a fully describable 
one-dimensional flow of point particles connected by ideal 
springs. Particle motion is constrained before reaching an 
inlet and after passing an outlet, and the system is shown 
to exhibit chaotic dynamics when particles are driven si- 
nusoidally after crossing the outlet. The outlet driving 
force continually adds energy to the system. No drag 
force is present, but energy is dissipated when particles 
exit at any point away from their equilibrium positions. 
The model can be reduced to a one-dimensional map that 
produces chaotic dynamics, showing that chaos can occur 
in flows at the boundary between simple constrained mo- 
tion, even when particles in the flow interact with linear 
forces. 

This work was supported by the National Science 
Foundation Grant No. NSF PHY 01-40179, NSF DMS 
03-25939 ITR, and NSF DGE 03-38215. 



Electronic address: gerig@uiuc.edu 

Electronic address: a-hubler@uiuc.edu 

K. R. Sreenivasan, Rev. Mod. Phys. 71, S383 (1999). 

D. Helbing, Rev. Mod. Phys. 73, 1067 (2001). 

D. Helbing, I. Farkas, and T. Vicsek, Nature 407, 487 

(2000). 

H. M. Jaeger, S. R. Nagel, and R. P. Behringer, Rev. 
Mod. Phys. 68, 1259 (1996). 

J. P. Gollub and J. S. Langer, Rev. Mod. Phys. 71, S396 



4 



(1999). 

[6] H. Aref and J. P. Gollub, in Research Trends in Fluid 
Dynamics, Report From the United States National Com- 
mittee on Theoretical and Applied Mechanics, edited by 
J. L. Lumley (American Institute of Physics, Woodbury, 
NY, 1996), p. 15. 

[7] S. Wolfram, J. Stat. Phys. 45, 471 (1986). 

[8] U. Frisch, B. Hasslacher, and Y. Pomeau, Phys. Rev. 
Lett. 56, 1505 (1986). 

[9] A. Brandstiiter, J. Swift, H. L. Swinney, A. Wolf, J. D. 
Farmer, E. Jen, and J. P. Crutchfield, Phys. Rev. Lett. 
51, 1442 (1983). 
[10] P. Holmes, J. L. Lumley, and G. Berkooz, Turbulence, 



Coherent Structures, Dynamical Systems, and Symmetry 
(Cambridge University Press, Cambridge, Great Britian, 
1996). 

[11] M. Peyrard and I. Daumont, Europhys. Lett. 59, 834 
(2002). 

[12] R. Z. Sagdeev, D. A. Usikov, and G. M. Zaslavsky, Non- 
linear Physics: From the Pendulum to Turbulence and 
Chaos (Harwood Academic Publishers, Chur, Switzer- 
land, 1988). 

[13] O. M. Braun and Y. S. Kivshar, Phys. Rep. 306, 1 (1998). 
[14] M. J. Feigenbaum, J. Stat. Phys. 19, 25 (1978). 



